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ABSTRACT 

We present a novel approach to derive the age of very young star clusters, by using the Turn-On 
(TOn). The TOn is the point in the color- magnitude diagram (CMD) where the pre-main sequence 
(PMS) joins the main sequence (MS). In the MS luminosity function (LF) of the cluster, the TOn 
is identified as a peak followed by a dip. We propose that by combining the CMD analysis with 
the monitoring of the spatial distribution of MS stars it is possible to reliably identify the TOn in 
cxtragalactic star forming regions. Compared to alternative methods, this technique is complementary 
to the turn-off dating and avoids the systematic biases affecting the PMS phase. We describe the 
method and its uncertainties, and apply it to the star forming region NGC346, which has been 
extensively imaged with the Hubble Space Telescope (HST). This study extends the LF approach in 
crowded extragalactic regions and opens the way for future studies with HST/WFC3, JWST and from 
the ground with adaptive optics. 

Subject headings: Magellanic Clouds — stars: formation — stars: pre-main sequence — galaxies: star 
clusters 



1. INTRODUCTION 

The MS turn-off is the most reliable feature for age- 
dating a star cluster. Nevertheless for very young clus- 
ters, with ages of tenths to few Myr, the identification of 
the turn-off is usually hampered by the paucity of mas- 
sive stars. 

We propose to take a different point of view, fo- 
cusing on the TOn, the CMD locus where the PMS 
joins the MS. Although the importance of the TOn 
has been already emphasized in s everal papers (e.g . 
Staufferl [19801 IBelikov et all [l99l IBaume et all [20031 
Navlor et al.H2009T) . its application to dating extragalac- 
tic star forming regions is a new proposition. 

In analogy with the turn-off, the TOn properties are 
directly related to the age of the stellar population, but 
with evolutionary times much shorter than the corre- 
sponding MS times. In fact, from simple stellar evolu- 
tion arguments, the age of a cluster is equal to the time 
spent in the PMS phase by its most massive star still in 
the PMS phase. By definition, this star is at the TOn. 
Hence, when the intrinsic luminosity of the TOn is de- 
tected, it is straightforward to associate it to the age of 
the cluster. 
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In the first part of this letter we describe how the in- 
trinsic properties of the TOn can be used as a clock. In 
the second part we present a new method to apply TOn 
related properties to date extragalactic systems. Finally 
we apply this method to the largest extragalactic star 
forming region, NGC346, in the Small Magellanic Cloud 
(SMC). 

2. THE PMS TURN-ON 

The potential strength of the TOn is apparent from 
the morphology of isochrones taking both the PMS and 
MS phases into account. FigQTa) shows 5 isochrones 
with metallic ity Z = 0.004, obtai ned combining the Pisa 
PMS tracks (ICignoni et al.ll2q09R w ith the Padua evolu- 
tionary tracks ([Fagotto et al.l[l994[ ) to cover the entire 
mass range 0.45 — 120 M Q . For all ages younger than 15 
Myr the isochrone portion just before the zero age MS 
(ZAMS) has a hook and then is significantly flatter than 
the ZAMS. The TOn is at the vertex of the hook, quite 
easy to recognize. 

To test if/how the TOn can be useful as a cosmo- 
chronomcter, we simulated synthetic simple stellar popu- 
lations (SSPs), using the tracks quoted above. As an ex- 
ample, we describe the case with metallicity Z = 0.004, 
burst duration 1 Myr, Salpeter's Initial Mass Function 
(IMF), and no binaries. To guarantee statistically signif- 
icant tests, extensive theoretical simulations have been 
performed with a number of synthetic stars up to 1 x 10 6 . 

Panels (b)-(f) in Fig [1] show the luminosity functions 
(LFs), binned in 0.2 magnitude bins and scaled to the to- 
tal number of stars, of synthetic SSPs (open histograms) 
corresponding to the isochrones of panel (a). A refer- 
ence synthetic zero age population, that is artificially 
built without PMS stars (grey filled histogram), is also 
shown. In the LF without PMS, the only valuable fea- 
ture is a mild peak around My « 2, consequence of an 

9 i.e. stars are born directly on the ZAMS. 
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Figure 1. Panel (a): combined PMS and MS isochrones for the labeled ages. The TOn is at the discontinuity where the PMS joins the 
MS (e.g. at (V — I)o = —0.15 and My = 0.85 for the 2 Myr case). The numbers on the left give the mass of the TOn star of each isochronc. 
Panels (b) to (f): normalized LFs for synthetic SSPs of the labeled age (see text for details). 



inflection in the derivative of the mass- My relation in MS 
stellar models around m = 2 M . By contrast, when the 
evolution starts from the PMS phase, the corresponding 
LFs develop an additional strong peak followed by a dip. 
Peak and dip reflect respectively the steep dependence of 
st ellar mass on My near the TOn (see also the discussion 
in lPiskunov et af1l2004D and the following flattening be- 
low the TOn (caused by the short evolutionary timescale 
of the PMS phase compared to the MS). After the dip, 
the shape of the LF mimics the IMF. 

The LFs in Fig |T] indicate the importance of both fea- 
tures, peak and dip, to infer the cluster age: the older 
the age, the fainter is the LF TOn peak and the corre- 
sponding dip. On the other hand, for the explored range 
of ages the magnitude of the MS peak is fairly constant 
(My « 2), being it locked to the (much longer) MS evo- 
lutionary times. Through a polynomial fit to the models, 
theory provides a useful relation between age (r) of the 



SSP and the magnitude My of the TOn in the range 
2-100 Myr: 



r{Myr) = E 7 =o,5%' x ( M v) j 



In turn: 



Mv 
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Once the bin width of the LF is chosen, Eq. (QJ (whose 
coefficients are given in Table [T]) gives also the intrinsic 
uncertainty on the age. Assuming 0.2 mag wide bins, 
the minimum uncertainty is about 0.6 Myr at 3 Myr, 1.3 
Myr at 20 Myr, 2.5 Myr at 30 Myr and 6 Myr at 50 Myr. 
The TOn formula is deliberately limited to populations 
older than 2 Myr, since for younger ages th e current PMS 
mode ls are still extremely uncertain (see iBaraffe et al.l 
120021 for a discussion). 

Although attractive, the TOn dating method should be 
used with caution. When we compare the observed and 
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Table 1 

Coefficients for Eq. 1 and Eq. 2. 
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theoretical LFs it is important to evaluate the following 
uncertainties: 

Incompleteness and photometric errors. To test real 
conditions, the synthetic SSPs have been degraded as- 
suming reddening and distance of the SMC, E(B — V) — 
0.08 and (m—M)o = 18.9, and photome t ric er rors and in- 
completeness as derived bv lSabbi et al.l (|2007|) from HST 
images of NGC346. Fig[2ja) shows the degraded LF for 
populations of 5, 30, 40, 50 Myr: at the distance of SMC 
the excellent quality of the HST/ ACS photometry guar- 
antees perfect detectability of the peak/dip feature up 
to about 50 Myr. Beyond this age the TOn becomes too 
faint to be detected. 

Poisson fluctuations. In order to check how many stars 
are necessary to safely identify the TOn, we progressively 
reduced the number of synthetic stars belonging to the 
5 Myr population, until the TOn peak was hidden by 
the Poisson fluctuations. This experiment indicates that 
about 50 stars brighter than My « 5 (corresponding to a 
cluster total mass of « 500 M Q ) are sufficient to identify 
the TOn with a significance of 2 a. 

IMF, binaries, star formation duration. In Figj2jb) 
the synthetic 5 Myr SSP has been modeled changing 
the binary fraction (50% of simulated stars have now 
an unresolved companion, randomly extracted from the 
same IMF as the primary) and the IMF exponent (1.5- 
3). None of these changes alters significantly the TOn 
position: the case with binaries is almost identical to the 
reference case, while the adoption of different IMFs mod- 
ifies only the shape of the LF before and after the TOn 
peak. 

In Fig. [2tc) we have simulated a 10 Myr old cluster 
with a st ar formation activity th at lasts 1, 3 and 5 Myr 
(see e.g. iPalla fc Stahleri l2000f ). As expected, with a 
prolonged SF activity the LF peak grows brighter and 
broader. In this case, the peak magnitude provides in- 
formation on the average age of the population. 

Reddening. Typical of star forming regions, the pres- 
ence of highly obscuring material (foreground as well 
as local) can significantly dim and blur a TOn. Thus, 
reliable reddening estimates are fundamental to obtain 
unbiased ages. On the other hand, there are several re- 
gions with affordable foreground and intrinsic extinction. 
NGC346 is one of them with foreground E(B-V) ~ 0.08 
and internal < 0.1 mag (as deduced from the upper main 
sequence), which contributes to the final uncertainty by 

1-2 Myr. 

3. IDENTIFICATION OF THE TON IN EXTRAGALACTIC 
STAR FORMING REGIONS: THE CASE OF NGC346 

Our goal is to study how star formation develops 
in extragalactic star forming regions. In the CMD of 
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Figure 2. Panel (a): Effect of incompleteness and photometric 
errors on the LFs for the labeled ages. All LFs are normalized 
to the total number of stars. Panel (b): Effect of different IMFs 
and binary fractions on the 5 Myr LF. The shadowed histogram 
is the same as the reference 5 Myr case of panel (a). Panel (c): 
normalized LFs for a 10 Myr cluster with the labeled duration of 
prolonged star formation activity. 
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Figure 3. The top-left panel shows the CMD selection (see text) of bona-fide MS stars (black dots). Each of the other panels shows the 
spatial distribution of these MS stars for the labeled range of magnitudes. 



these regions, the MS is often contaminated by young 
fore/background stars and only in a few cases member- 
ship information is available for a safe decontamination. 
The presence of the field MS partially fills the LF dip, 
and sensibly lowers the significance of the TOn in the 
LF. 

To get around the problem of contamination we pro- 
pose to combine the fact that, by definition, no sub- 
cluster member on the MS can be fainter than the TOn, 
together with a careful analysis of the spatial distribution 
of the stars in the region. The procedure we suggest has 
3 main steps: 1) Selection of bona-fide MS stars (both 
members and non-members of the clusters) of all magni- 
tudes. In order to take into account photometric errors 
and reddening, we consider all the stars bluer than a ridge 
line appropriately redder than the theoretical ZAMS. 2) 
Division of the selected MS stars into bins of progres- 
sively fainter magnitude. For each bin of a given mag- 
nitude, the spatial distribution of stars is examined to 
assess whether the sub-clusters are still visible. When a 
sub-cluster is clearly identified up to a given magnitude, 
but it disappears in the fainter maps, we identify the 
magnitude of its TOn. This is because we have reached 
the magnitude where the sub-cluster members have not 
yet reached the MS, and MS stars fainter than this limit 
do not belong to that sub-cluster. In order to evaluate 
the statistical significance of the sub-cluster disappear- 
ance from the maps, for each of the three sub-clusters 
we built a stellar density radial profile and we evaluated 
the contribution from different magnitude bins. For field 
stars we expect a flat stellar density profile, while a de- 
crease in stellar density is the typical signature of a sub- 
cluster. The range of magnitudes where the transition 
between the two profiles occurs identifies the apparent 
magnitude of the TOn. 3) Age determination of the sub- 
clusters. After applying reddening and distance modulus 
corrections, we use Eq. ([T]) to translate the sub-cluster 
TOn magnitude into an age. The final uncertainty is 
fixed by the bin width, which is, in turn, chosen to ob- 



tain an acceptable number of counts per bin. 

To test the strength of our method, we applied it to 
the star forming region NGC346 in the SMC, whose im- 
ages acquired by the HST Advanced Camera for Survey 
have revealed a wealth of young sub-clusters containing 
several PMSs (ISabbi et al J 120071) . This region provides 
excellent conditions to test our method because its recent 
stron g activity supplie s an outstanding sample of PMS 
stars (Nota et al. 2006). The complexity of its structure, 
with several non coeval sub-clusters, requires a method 
able to trace the temporal sequence of events leading to 
the present configuration. In this letter we focus on three 
of the ri chest sub-clu s ters, SC-1, SC-13 and SC-16, as de- 
fined bv ISabbi et~aT1 (|2007j ). Their location in the region 
can be seen in their Fig. 8. 

Fig|3] summarizes our analysis of NGC346. To select 
only MS stars we used a ridge line 0.05 mag redder than 
the theoretical Z = 0.004 ZAMS. Wc considered all the 
stars bluer than this line as bona-fide MS stars, and di- 
vided them in six bins of different size. Bona-fide MS 
stars are drawn in black in the top-left panel of FigJ3J 
where all the other stars are marked in grey. The other 
five panels show the spatial distribution of the bona-fide 
MS stars for progressively fainter bins of magnitude. An- 
alyzing these maps, we find that the central sub-cluster 
SC-1 is well visible down to V = 21, while no obvious 
spatial structure resembling SC-1 is observable in the 
fainter maps. The sub-cluster SC-16, on the contrary, is 
still clearly visible down to V — 23. Beyond this limit, 
also SC-16 vanishes. The small sub-cluster SC-13 is rec- 
ognizable at least to V = 21. Notice that the stellar 
agglomeration appearing at V = 22 corresponds to the 
4-5 Gyr turn-off s tars of the older cluster BS90 (see e.g. 
ISabbi et al J 12001 in the foreground of NGC346. 

In order to test the significance of the sub-clusters 
we examined the radial density profile of the three sub- 
clusters as a function of magnitude (Fig. [4] left pan- 
els). The profiles are calculated using annuli of equal 
area centered on the highest density peak. The inner 
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radius is fixed at 150 pixel. To exclude any complete- 
ness issue in these crowded regions, we also show in the 
middle panels of Figure [4] the completeness factors ob- 
tai ned from the exte nsive artificial star tests performed 
bv lSabbi et al.1 (j2007j ). for all stars within 100 pixels from 
the centers of SC-16, SC-1 and SC-13. To further test 
the final ages, theoretical isochrones are over-imposed on 
the CMDs (right panels of FiggJ). 

3.1. SC-16 

The first row in Fig0]shows the MS radial profiles (left 
panel), the completeness curves (middle panel) and the 
CMD (right panel) for SC-16. The compact morphology 
of this sub-cluster is evident: the radial profile of MS 
stars brighter than V = 23 rapidly drops with distance 
from the center, while the distribution of less luminous 
MS stars is constant. This confirms that the TOn magni- 
tude is between 22 and 23, thereby constraining the first 
generation of stars to have formed between 12.5 and 18 
Myr ago. The upper right panel of Fig|4] shows the cor- 
responding isochrones superimposed to the SC-16 CMD 
(all stars within 150 pixel from the center): it is reassur- 
ing to see that below the 18 Myr TOn, the sub-cluster 
stars move away from the MS, confirming that the MS 
deficit in the map is a genuine evolutionary effect. Con- 
cerning the completeness, although the sub-cluster region 
loses stars faster than the field, for magnitudes V < 23 
this effect is not severe (compared to the field, less than 
25% of stars are lost). By comparing our age estimate 
to indepe ndent determ i nation s we find good agreement 
both with ISabbi et all (120071) (15 ± 2.5 Myr) and with 
IHennekemper et al.1 (|2008l ) (5-15 Myr). 

3.2. SC-1 

The case of SC-1 (second row in Fig. HJ) is quite dif- 
ferent. Here the extended spatial morphology produces 
a gradual decrease in the MS radial profile rather than 
an abrupt drop at V 21. Radial profiles for fainter 
magnitudes are quite flat. Assuming a conservative bin 
uncertainty of 1.5 mag on the TOn magnitude, the star 
formation onset is expected between 3.5 and 6.5 Myr ago , 
in agreement with the estimate by Sabbi et al.l (|2007T ) 
(3±1 Myr). Both the 3.5 and the 6.5 isochrones are over- 
laid on the SC-1 CMD. As suggested by radial profiles at 
V = 21 there is a clear signature of a PMS TOn, with a 
large sample of stars still in the PMS phase. Our 6.5 Myr 
isochrone fits the upper MS morphology, including the 
TOn, very well and fits the PMS blue envelope as well. 
The few MS stars below the TOn are likely contamina- 
tion of the SMC field and/or stars of the old cluster BS90 
(the presence of 3 clump stars endorses this hypothesis). 
The absence of any star on the MS around V ~ 21.5 and 
the morphology of the CMD for the vast majority of the 
stars, which seems to follow the isochrone with PMS in- 
cluded, further supports our conclusion. The measured 
completeness (middle panel) is better than 75% and can- 
not account for such a reduction. 

3.3. SC-13 

The third row of Fig. @] shows the sub-cluster SC- 
13. As for SC-1, the MS density profiles of SC-13 arc 
declining down to V > 21. However, due to the poor 
statistics, we can only state that the cluster formation 



started sometimes in the last 6 Myr. There is a fur- 
ther clue suggesting that SC-13 is actually younger than 
6 Myr: a close inspection of the CMD (right panel) 
reveals that most of probable intermediate-mass PMS 
stars are aligned with the 3 Myr isochrone. In agree- 
ment with this finding, most of the low-mass PMS stars 
are redder th an the 3 Myr iso chrone. This age is 
also found by ISabbi et al l (|2007t ) (3 ± 1 Myr) and by 
IHennekemper et al.l (|2008f ) (0.5-2.5 Myr). 



4. DISCUSSION AND CONCLUSIONS 

The evolution of young star forming regions in their 
earliest stages is still quite unknown. For very young 
clusters, with ages of tenths to few Myr, the identifica- 
tion of the turn-off is usually hampered by the paucity 
of massive stars. However when a cluster is sufficiently 
young to harbor PMS stars, the luminosity of its TOn 
provides a robust indication of its age. Furthermore, 
while the PMS evolutionary models still suffer of several 
limitations, the simple comparison of TOn luminosities 
is a reliable measure of relative ages. 

In the LF the TOn is a narrow peak followed by a dip. 
The TOn is very sensitive to age: for a <~ 30 Myr old 
stellar population the luminosity of the TOn changes by 
~ 0.08 mag/Myr, but at ~ 20 Myr it already changes by 
~ 0.15 mag/Myr, and by ~ 0.33 mag/Myr at 3 Myr. 

To guarantee a safe TOn identification, it is impor- 
tant to select targets with low or well known reddening. 
In the visual, this currently excludes Galactic clusters 
like Westerlund 1-2, NGC3603, Arches and Quintuplet, 
but good targets can be found in the solar vicinity and 
in some nearby irregulars galaxies like the Magellanic 
Clouds and IC1613. 

In this letter we have presented a new method to in- 
vestigate how star formation develops in complex extra- 
galactic young clusters such as NGC346 in the SMC. 

Our method combines the analysis of the star clusters 
stellar densities profiles with the notion that in a cluster 
no star on the MS can be fainter than the cluster TOn. 
This approach has the advantage of strongly reducing 
the uncertainties introduced by the contamination to the 
CMD by young stars that do not belong to the cluster, 
ultimately affecting the TOn detectability in both CMD 
and LF. 

Clearly, there are limitations to the applicability of the 
method. The bona-fidc MS selection may include intrud- 
ers like PMS and stars evolved off the MS, reducing the 
TOn visibility. The issue is particularly thorny at ages 
larger than about 30 Myr, because the PMS isochrones 
are close to the ZAMS (see Figure QJa) ) ■ The method 
relies on the assumption that the formation sites are 
agglomerations of stars: if for some reason the new- 
born stars formed in isolation or, drifting apart, were 
too diluted to be detected as aggregates, the method 
would not be a pplicable. In this respect, the study by 
Pfalzncr (2009) , exploiting the density-radius relation to 
date nearby clusters younger than 20 Myr is interesting 
and reassuring. 

We applied the TOn method to NGC346 in SMC, and 
we found that the onset of the star formation in the sub- 
cluster SC-1 was between 3.5 — 6.5 Myr, in the sub-cluster 
SC-13 3 Myr ago or less and in SC-16 between 12.5 — 18 
Myr ago, in good agreement with the results from the 
literature. 
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Figure 4. Left-hand panels: radial distribution of MS stars near the centers of SC-16 (top), SC-1 (middle) and SC-13 (bottom). Different 
curves correspond to the indicated ranges of magnitude. The abscissa is in units of n, where n is related to the distance d from the 
sub-cluster center by d = \fn X 150 pixels. Middle panels: Completeness for the labeled sub-clusters (dashed line) and field stars (solid 
line). Right-hand panels: CMDs for all stars within 150 pixels from the center of SC-1 (top panel), SC-16 (middle panel) and SC-13 
(bottom panel). Red dots indicate sub-cluster stars while black dots stand for the entire data sample. Isochrones of the indicated ages arc 
also shown. 



Having established the effectiveness of our method, the 
next steps will be to: 1) incorporate near-infrared pho- 
tometry and use it to study extinguished regions where 
measurements of the upper MS are difficult to carry out; 
2) undertake comparative studies of the duration and 
spatial patterns of SF in Magellanic Cloud regions rang- 
ing from the giant 30 Dor complex to the comparatively 
isolated NGC602 cluster. 
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